Slowly Rotating Black Holes in Dynamical Chern-Simons Gravity: 
Deformation Quadratic in the Spin 
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We derive a stationary and axisymmetric black hole solution to quadratic order in the spin angular 
momentum. The previously found, linear-in-spin terms modify the odd-parity sector of the metric, 
while the new corrections appear in the even-parity sector. These corrections modify the quadrupole 
moment, as well as the (coordinate-dependent) location of the event horizon and the ergoregion. 
Although the linear-in-spin metric is of Petrov type D, the quadratic order terms render it of type I. 
The metric does not possess a second-order Killing tensor or a Carter-like constant. The new metric 
does not possess closed timelike curves or spacetime regions that violate causality outside of the 
event horizon. The new, even-parity modifications to the Kerr metric decay less rapidly at spatial 
infinity than the leading-order in spin, odd-parity ones, and thus, the former are more important 
when considering black holes that are rotating moderately fast. We calculate the modifications 
to the Hamiltonian, binding energy and Kepler's third law. These modifications are crucial for 
the construction of gravitational wave templates for black hole binaries, which will enter at second 
post-Newtonian order, just like dissipative modifications found previously. 

PACS numbers: 04.30.-w,04.50.Kd,04.25.-g,04.25.Nx 



o 
cr 



(N 
> 
O 
m 

O 

(N 



X 



I. INTRODUCTION 

Although General Relativity (GR) has passed all Solar 
System and binary pulsar tests it remains yet to be 
verified in the non-linear, dynamical strong- field/strong- 
curvature regime. One of the best astrophysical environ- 
ments to perform such tests is near black holes (BHs). 
Mathematical theorems in GR guarantee that the exte- 
rior gravitational field of a vacuum, stationary and ax- 
isymmetric BH is given by the Kerr metric [2t2| The 
multipolc moments of an uncharged BH in GR are there- 
fore completely characterized by its mass (mass monopole 
moment) and spin angular momentum (current dipole 
moment). If GR is modified, however, BHs need not be 
represented by the Kerr solution. For example, a sta- 
tionary BH solution that differs from the Kerr metric 
has recently been found in quadratic gravity to linear- 
order in a slow-rotation approximation [8l4ll|. Future 
observations of electromagnetic radiation from accretion 
disks fl2l4h3 | and gravitational radiation from extreme 
mass-ratio inspirals (l5rl2~j| will allow us to probe the 
spacetime structure and test the Kerr hypothesis: that 
the massive compact objects at the center of most galax- 
ies are Kerr BHs. 

A well-motivated theory that has recently received 
much attention is Chern-Simons (CS) modified grav- 
ity 0, [23| . The CS action modifies the Einstein- Hilbcrt 
one by adding a kinetic scalar field term and a potential 
composed of the product of this field with the Pontrya- 
gin density (the contraction of the Riemann tensor with 
its dual) . Such a potential is necessary to cancel anoma- 
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lies in heterotic, superstring theory |24| and it also arises 
naturally in loop quantum gravity [25l427j and in effective 
field theories of inflation [28| . This theory is also related 
to 3D topologically massive gauge theories [29| . Two for- 
mulations of CS gravity exist: a non-dynamical one and 
a dynamical one. In the former, the scalar field kinetic 
term is absent from the action, and thus, this field must 
be prescribed a priori, leading to an extra constraint on 
the space of solution, i.e. the vanishing of the Pontryagin 
density [3(|. In the latter, the scalar field kinetic term is 
kept in the action and the scalar is controlled by a wave 
equation sourced by the Pontryagin density. Therefore, 
the dynamical theory preserves diffeomorphism invari- 
ance and the strong equivalence principle, although it vi- 
olates GR's Birkhoff theorem [3(J, [HJ and the effacement 
principle (32[. 

Dynamical CS gravity is more appealing from a the- 
oretical standpoint but it remains relatively unexplored 
due to its mathematical complexity. The first study of 
BH solutions in dynamical CS gravity was carried out by 
Yunes and Pretorius 0, and later by Konno et al. 
In these studies, a slow-rotation approximation was used 
to find the leading-order in spin corrections to the Kerr 
metric, which enter in the gravitomagnetic sector and 
modify frame-dragging. A rather weak but robust con- 
straint on dynamical CS gravity was placed [33j j with 
Gravity Probe B [34[ (similar constraints can be ob- 
tained from LAGEOS and LAGEOS 2 [pi through mea- 
surements of the Lense-Thirring effect |36[.) Dynamical 
CS gravity can also be constrained by table-top experi- 
ments that confirm Newton's law to length scales above 
O(10^m) [13] ■ Interestingly, table-top experiments place 
similar constraint to those obtained from Solar System 
observations. 

A much stronger constraint can be placed with future 
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gravitational wave observations of compact binary coales- 
cences. This necessitates a calculation of both dissipative 
corrections to the waveform, i.e. modifications to the rate 
of change of the energy, angular momentum and Carter 
constant, and conservative corrections, i.e. modifications 
to the binding energy and Kepler's Third Law. In [3^ |. 
we calculated the former and found that they introduce 
modifications to the waveform at 2nd post-Newtonian 
(PN) order in terms of its frequency (or velocity) de- 
pendence, although its amplitude is of course suppressed 
relative to GR by the CS coupling constant. In order 
to obtain the full 2PN waveform, we must first compute 
the 2PN conservative corrections, which require a calcu- 
lation of (i) a stationary, axisymmetric BH solution at 
quadratic order in spin, and (ii) the dipole-dipole scalar 
force acting on the binary components. In this paper, we 
concentrate on the former. 

We will search for a slowly-rotating BH solution in dy- 
namical CS gravity at quadratic order in spin with two 
approximations: slow- rotation and small coupling Q. 
The former assumes that the BH spin angular momen- 
tum is much smaller than its mass squared. The latter 
postulates that the deformation away from GR is small, 
which corresponds to a dynamical CS dimensionless cou- 
pling constant [defined later in Eq. ([I])] much smaller than 
unity. This is a reasonable approximation, given that 
GR has already passed stringent tests, albeit in the weak 
field. Moreover, the dynamical CS action is an effec- 
tive theory that derives from a leading-order truncation 
in the dynamical CS coupling parameter of a more fun- 
damental one. Therefore, the action and its associated 
modified field equations are only valid to leading-order 
in the dynamical CS coupling parameter. If one did not 
use the small-coupling approximation to itcratively solve 
the modified field equations, third time derivatives could 
arise in the field equations, probably signaling the pres- 
ence of ghost modes. Presumably, these ghost modes 
might be eliminated by an appropriate UV completion. 

The method we employ to find a slowly-rotating solu- 
tion at second order in the spin is rather novel in this 
context. We follow closely BH perturbation theory tech- 
niques [H, Hi! : we treat the second-order in spin correc- 
tion to the dynamical CS metric as a perturbation away 
from the leading-order in spin solution found by Yuncs 
and Pretorius [g. We then derive master equations gov- 
erning this perturbation and decompose the solution in 
tensor spherical harmonics, which decouples the system 
into a linear ordinary differential one. The solution to 
this system is then mathematically straightforward, al- 
though we do verify that the full metric satisfies the dy- 
namical CS field equations explicitly with symbolic ma- 
nipulation software. As expected, the quadratic-order in 
spin corrections modify only the even-parity sector of the 
metric. 

The properties of this solution are also quite interest- 
ing. Since we obtain corrections to quadratic order in the 
spin, we are able to compute the gauge-dependent shift in 
the location of the event horizon and the ergosphere, as 



well as the deformation in the Kerr quadrupolc moment. 
We also find that the new metric retains its Lorentz signa- 
ture and closed timelike curves do not exist outside of the 
event horizon. Therefore, such a spacetime is perfectly 
suitable to the study of photon trajectories when consid- 
ering the shadow of BH accretion disks. We also consider 
test-particles motion in this new background. We obtain 
the corrections to the relation among the binding energy, 
the (z-component of) angular momentum, and the orbital 
frequency, and the frequency at the innermost stable cir- 
cular orbit. 

Perhaps one of the most interesting, albeit not sur- 
prising properties of the new solution is that, although 
the linear order in spin metric @ remains of Petrov type 
D [4(|, the quadratic order terms calculated here ren- 
der the metric of type I. Obviously, this is drastically 
different from the Kerr metric, which remains of Petrov 
type D to all orders in a small spin expansion. This 
is because if a metric is not of Petrov type Dat some 
order in spin, this violation cannot be compensated for 
by higher order spin contributions, as long as the small 
spin approximation holds. The new BH solution does 
not possess a non-trivial second-rank Killing tensor, and 
thus, there are no naive extensions of the Carter-like con- 
stant and a separable structure is not admissible fill . |42| . 
The latter statement means that there does not exist any 
spacetime coordinate transformation that leads to sepa- 
rable Hamilton- Jacobi equations. One cannot, however, 
rule out a canonical transformation that does render the 
equations separable. This is related to the fact that there 
might exist a higher-rank Killing tensor that kneads to a 
fourth constant of the motion. Given this, it is not clear 
whether geodesies in this new background will be chaotic 
or not. Nonetheless, upon orbit-averaging, a modified 
Carter-like constant reappears, which might indicate that 
orbits are regular, except for resonant ones. 

The remainder of this paper presents further details 
and it is organized as follows. Section |TT] presents the 
basic equations of dynamical CS gravity. Section |HT] be- 
gins by describing the approximation used to find BH 
solutions in this theoryand then continues to describe 
the solution found in [8| and @ and the new solution 
found in this paper. Section IIVI investigates the basic 
properties of this new BH solution, such as the (gauge- 
dependent) location of the event horizon and ergosphere. 
SectionfV] discusses geodesies in this new background and 
the Petrov type of the new solution. Section I VII sum- 
marizes and describe several possible avenues for future 
research. 

All throughout, we use the following conventions, fol- 
lowing mostly Misner, Thorne and Wheeler (43|. We 
use the Greek letters (a, /?,•■■ ) to denote spacetime in- 
dices. The metric is denoted g^ v and it has signature 
(— , +, +, +). We use geometric units, with G = c = I. 
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II. DYNAMICAL CHERN-SIMONS GRAVITY 

In this section, we introduce the basic equations of 
dynamical CS gravity and establish notation. The action 
is defined by [23( 



5* = 



;[V^V^ + 2V(tf)]+£ mat }. (1) 
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Here, n g = (167r) _1 , g denotes the determinant of the 
metric g^ and R^ v scr is the Riemann tensor. *R^ u P a is 
the dual of the Riemann tensor, defined by [23| 



(2) 



where ^ va P is the Levi-Civita tensor. t9 is a scalar field, 
while a and /3 are coupling constants. V^tf) is an addi- 
tional potential and £ ma t denotes the matter Lagrangian 
density. Following ||, we take $ and /3 to be dimension- 
less and a to have dimensions of (length) 2 . For conve- 
nience, we define a dimensionless parameter 



a 



KgfiM 4 



(3) 



where M is the typical mass of the system (or strictly 
speaking, the curvature length scale). 

The field equations in this theory are given by [23| 

^ + fc,„ = -^(T™ t +T; i ,) ) (4) 

Kg ZKg 

where G^ IV is the Einstein tensor and T™ at is the matter 
stress-energy tensor. The C-tensor and the stress-energy 
tensor for the scalar field arc defined by 



Equivalently, Eq. ([4]) can be rewritten as 



2k„ 1 



T" 9 ) 



(5) 

(0) 
(7) 



where we introduced the trace-reversed stress-energy ten- 
sors 



1 



rriT. 



rpti 



/3(V M i?)(V I/1 ?) 



(8) 
(9) 



and used the fact that the C-tensor is traceless. The 
evolution equation of the scalar field is given by 



4/3 



±Xj VilQ(7 Jt 



dV_ 



(10) 



In 4D, low-energy effective string theories, the (axion) 
scalar field has a shift symmetry, i.e. the equations of 
motion are invariant under the symmetry transformation 
i9 — > $ + c, with c a constant, which disallows any mass 
terms in the action. If one forgets about shift-symmetry 
and insists on a mass term, then the dynamics of the 
scalar field would freeze and their would not be much 
scalar-field propagation. In dynamical CS gravity, how- 
ever, such mass terms are not allowed. For simplicity, in 
this paper we set V($) = throughout. 

Dynamical CS gravity should be thought of as an effec- 
tive theory, and as such, it possess a cut-off scale outside 
which its action should be modified through the inclusion 
of higher-order curvature terms. This cut-off scale can be 
determined by estimating the order of magnitude of loop 
corrections to the second term in Eq. (fT]) due to rt-point 
interactions. Denote the additional number of vertices 
and scalar and graviton propagators relative to tree level 
diagrams as V, P s and P g , one has 



P, 



V 
2"' 



Pa 



1) 



V 



(11) 



One immediately sees that loop corrections arc sup- 
pressed by a factor of a v M^ nS>V K nV , where M v \ is the 
Planck mass and A is the energy scale introduced such 
that the suppression factor becomes dimensionless. This 
factor becomes of order unity when A takes the critical 
value 



A c = My' n oV n 



(12) 



which corresponds to the cutoff energy scale above which 
one cannot treat dynamical CS gravity as an effective 
theory. 

Given this cut-off scale, one can estimate the value of 
a above which the theory is not an effective one any- 
more. First, notice that for a fixed value of a, A c be- 
comes larger as n increases. Hence, n — 3 will lead to 
the most stringent constraint on a. From Eq. (|12p . the 
critical wavelength scale below which the strong coupling 
effect cannot be neglected is given by 



(13) 



where L p \ is the Planck length scale. When one takes A c 
to be of O(10um), thus saturating the table-top experi- 
ment bound 1371 . one finds that 



fa < O(10 8 km). 



(14) 



For values of a that satisfy the above inequality, dy- 
namical CS gravity can be treated as an effective the- 
ory and higher-order curvature terms in the action can 
be neglected. Notice that this inequality is of the same 
order as the constraint obtained from Solar System ex- 
periments [33]. Of course, one can have a value of a 
that satisfies this inequality, without necessarily having 
(<1. In this paper, however, we are interested in BHs 
of masses in the range (3,10 7 )M Q , for which the small 
coupling approximation requires yfo. <C 10 7 km, thus au- 
tomatically satisfying the inequality in Eq. (|14[) . 
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III. ROTATING BLACK HOLE SOLUTIONS 

In this section, we describe the two approximation 
schemes that we use to obtain the slowly rotating BH 
solution in dynamical CS gravity at quadratic order in 
spin. We then describe the slowly-rotating BH solution 
at linear-order in spin, found by Q and Q, and apply 
these approximations to find the second-order in spin so- 
lution. 



A. Approximation Schemes 



Let us now expand the scalar field d. From Eq. (|10l) . 
we see that the leading-order contribution to i? is propor- 
tional to a. Therefore, we can expand $ as 



x 'tf(i,D + x *0(2,i) + 0(a' X ' 3 ). 



(21) 



There is no i^ 0,1 ) term here because the Pontryagin den- 
sity vanishes when evaluated on spherically symmetric 
spacetimes. There is no 0(a' 2 ) term and we have also 
here neglected terms of 0(a' 3 ) since they do not affect 
the metric perturbation at 0(a' 2 ). 



Following we consider stationary and axisymmet- 
ric BH solutions in dynamical CS gravity with small- 
coupling (£ <C 1) and slow-rotation (x <C 1). The small- 
coupling approximation implies that we consider small 
CS deformations away from GR. The metric can then be 
expanded as 



V=.^+«' 2 ^+0(a' 4 ), 



(15) 



where a' is a bookkeeping parameter that labels the order 
of the small-coupling approximation, with gj$ oc a n . 
Notice that only terms with even powers in a' appear 
in the metric. Then, we expand each g$ and g 2 } in a 
slow-rotation expansion via 

*' 2 9$ = «' 2 9$ 2) +«' 2 x' 9 fr 2) +*V9% 2) 

+0(a' 2 x' 3 ) , (17) 

where x is another bookkeeping parameter that labels 
the order of the slow-rotation approximation. Notice that 

gflu oc x m a", where x = a/M is the dimensionless spin 
parameter. 

The quantities g^v , gfy and g^' ' can be obtained 
by expanding the Kerr solution in x ^ 1: whose line 
element in Boyer-Lindquist (BL) coordinates (t, r, 9, cj>) 
is 



B. BH solutions to £>(a' 2 x') 

Let us first concentrate on solutions at 0(a' 2 x' )- As 
already mentioned, the Pontryagin density vanishes for 
any spherically symmetric spacetime [30j. Thus, static, 
spherically symmetric BHs are still described by the 
Schwarzschild solution. This implies that g^v = for 



(0,2) 



0. 



all n, and in particular, 

To 0(a' 2 x') m metric, Yunes and Pretorius found that 
in BL-type coordinates (the coordinates where the GR 
part of the BH metric is identical to Kerr in BL coordi- 
nates) nm 



0(1,1) = 5ocob< 
8/3* r 2 



2 M 18 M' 



1+ + 
7 r 5 



-T ( 22 ) 



and the only non- vanishing term in gfy 2 ^ is 



(1,2) 

9t<p 



12 M 27 M 2 \ 

! + t— + 7?;— sin ' 23 

7 r 10 r 



with all other components set to zero. Therefore, the 
line-element to 0(a' 2 x') is given by 



ds 2 



ds 2 K 



5 ^r 



12 M 27 M 2 \ . 2ni 

1 + YV + iQ—) smddtd ^ 

(24) 



dsi 



2Mr\ , 9 4M ar sin 2 8 
1 dt 2 



+Y,dd 2 



1 i 2 

r + a 



dtdcj) + —dr 2 



2Ma 2 r sin 2 9 



sm 



(18) 



Notice that although the correction term does not diverge 
at the unperturbed Schwarzschild horizon, r = 2M, it 
does not vanish there, either. One can, however, resum 
the metric such that the correction indeed vanishes at the 
Schwarzschild horizon, as we will discuss in Sec. MI C 31 



where A and S are defined by 



A = r 2 - 2Mr + 
S = r 2 + a 2 cos 2 I 



(19) 
(20) 



Here, M is the mass of the BH and a = S/M with S 
denoting the magnitude of spin angular momentum of 
the BH. 



C. BH solutions at 0(a' 2 X ' 2 ) 
1. Scalar Field 



Since the right-hand side of Eq. (flOf is already pro- 
portional to a/P, we expand the Pontryagin density 
Ru tJ .pcr*R fnJf " J only up to 0(a' ). This means that we can 
substitute the Kerr solution in R vtipa * R^ upa and expand 
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it in powers of x' to find 

* , „„„ M 3 X cos Of 28 M 2 9 9 
Ru^pa R p = 288 ^ ( f - —— X 2 cos 2 i 



+o(x' 5 ) 



(25) 



The Pontryagin density is a parity odd quantity, and 
hence it can only depend on odd powers of x' ■ Since 
we are only interested in a BH solution to 0(a' 2 x' 2 ), 
we need not concern ourselves with the second term in 
Eq. (|2"5)) . Therefore, one finds that 



(2,1) =0 . 



(26) 



and we only have to consider tfC 1 ' 1 ) [see Eq. (|22p] to con- 
struct gfy 2 ^ ■ In fact, this shows that i^™' 1 ) = for all 
even n. 



2. Metric Tensor: Equations 

Consider an expansion of the metric of the form g^ v = 

9^v + where denotes a metric perturbation 
away from Schwarzschild solution. For us, this metric 
deformation contains both known terms, such as pure 
Kerr deformations of Schwarzschild and CS corrections 
at 0(a' 2 x'), as well as unknown terms, such as CS cor- 
rections at 0(a' 2 x' 2 ), namely 

Ku = X'9^ + X' 2 ^ + X'ol 2 ^ + x' 2 a' 2 g^ . 

(27) 

The Einstein tensor can then be expanded as 
GW = Gfl + GW [V] + Gf v [V. KA +0(h 3 )- (28) 

Here, the superscript in square brackets counts the num- 
ber of times appears. Obviously, the first term in 
Eq. (|28p vanishes because the Schwarzschild metric sat- 
isfies the vacuum Einstein equations. 

With this notation, we can split the 0(a' 2 x' 2 ) part of 

(2 2). 

the Einstein tensor G)iu into two contributions 



For convenience, we introduced the reduced C-tensor C'^ v 
and the reduced stress energy-momentum tensor of the 



scalar field T, 



■d /. 



ee (a/ Kg )C^ , T» ee (1/2^)2* 



(32) 



The components of S^J 2 * 1 can be calculated from 
Eqs. ([2"2")l. P3| and the Kerr metric. 

The rccastcd field equations carry a strong resemblance 
with the equations of BH perturbation theory [38l . l39j . 
The quantity on the left-hand side of Eq. ([3H)l can be in- 
terpreted as the Einstein tensor constructed from the un- 
known perturbation g^ in a Schwarzschild background 



(0,0) 



Since the source S 



(2,2) 



is an analytic function that 

can be computed exactly, g y2 J 2 ^ can be solved for using 
Schwarzschild BH perturbation theory tools. 

Following Refs. jH, H3, we first decompose the metric 



perturbation gfy 2 ^ and the source Sfxh> in tensor spher- 
ical harmonics. Since terms of 0(a' 2 x' 2 ) are parity even, 
we only consider the metric perturbation in the even- 
parity sector, which has 7 independent metric compo- 
nents. Imposing stationarity and axisymmetry reduces 
the number of independent components to 5. The latter 
condition allows us to consider the m = mode only in 
the spherical harmonic decomposition. Two gauge de- 
grees of freedom remain, which we fix by adopting the 
Zerilli gauge. One is then left with 3 independent de- 
grees of freedom, which allows us to parameterize the 
metric perturbation via 



? (2,2) 



,(2,2) 



+V2K eo (r)geo 



(33) 



and the source term via 



S (2,2) = [A$a$ + Ataato + B m b m 



n (2,2) = n [l] 



7 (2,2) 

J fits 



G 



„(1,0) (1,2) 



(29) 



(34) 



where the first term is constructed from the unknown 

(2 2) 

functions g^ and its derivatives, while the second term 
is a known function built from gfy°^ and gfy only. We 
can then rewrite the field equations at 0(a' 2 x' 2 ) as 



gw 



„(2,2) 



- C?( 2 > 2 ) 



(30) 



where we have defined the source term 



7(2.2) 



w (11/ 



(1,0) (1,2) 



pi (2,2) , rpl-Q (2,2) 



(31) 



.(0) 



1 — 2M/r and the five tensor spherical har- 
a i0i bio, 9 io and fm defined in Appendix [A] 



with f{r) 
monies a 

Notice that boldfaced quantities here refer to rank-2 co- 
variant tensors. The source term coefficients A^ \ Aiq, 

(s) 

Bio, G iQ and are non- vanishing only for £ = and 
I = 2, and we provide explicit expressions for them in 
Appendix [A] 

Using this decomposition, the field equations in 
Eq. (f3TJ|) are no longer coupled, partial differential equa- 
tions, but they rather become coupled ordinary differen- 
tial equations for (H^o, H2m, Km), namely [38l. l39j 



6 



f(r) 



, d 2 Ki, 
dr 2 



+ zf(r) 3 



5M\ 



dK 



CO 



1 

r J dr r 
1(1 + 1) 



fir) 



dH 



2/0 



dr 



Co) 



r^ " ' ' 2r 2 
r - M dK m , 1 dff />o , 1 



r 2 /(r) dr 



r dr 



r 2 f{r) 



f(r)(H 2m + K m ) = —A 



f(r)-^(H om - Kta) + + (H 2£0 - H m ) = / /(r) 5, 

dr r 2 r \ r J ^/JQ + l)/2 



- 



r + MdH om 1(1 + 1 



M\ dK eo ,d 2 H oeo 1 



/(r) 



dr 2 r y r J dr dr 2 r \ r J dr 



1 



M\ dff 



2m 



dr 



0£0 



2(0 



2r 2 
r 2 F eo 



(H eo — H 2 K)) = V^G^q , 



y/l(l + l)(l-l)(l + 2)/2 



(35) 
(36) 
(37) 

(38) 
(39) 



In Eqs. (35]), (36]) and (38]), ^ = or 2, while in Eqs. (37]) 
and (39]) £ = 2. Wc have checked that these equations are 
identical to those derived from the field equations directly 
through the use of symbolic manipulation software. 

3. Metric Tensor: Solution 

Before solving these equations, let us exhaust the re- 
maining gauge freedom in the i = mode. We already 
explained that for modes with I > 2, one is left with 
5 independent variables after imposing stationarity and 
axisymmetry. Mathematically, these variables are con- 
tained in the coefficients of the five spherical tensor har- 
monics 0!$ , ago, bio, 9eo an d feo- Wc eliminate two 
of them by imposing the Zerilli gauge, i.e. setting the 
coefficients of beo and fio to zero. The I = mode, 
however, is different because from the start it possesses 
only 3 independent variables, after imposing stationary 
and axisymmetry. One of them corresponds to a redef- 
inition of the spherical areal radius, which we eliminate 



by setting K o = 0. 

Let us now discuss how to solve the differential sys- 
tem in Eqs. (33 -(33). When wc substitute K 00 = in 
Eq. (|35j) with I = 0, we are left with a first-order or- 
dinary differential equation for H 2 oo- We can solve for 
H 200 and then use Eqs (J3BJ) and (3E} to find H 000 . With 
this, we can then obtain the I = 2 perturbative modes 
to find the (Hoeo, H 2 £o, Km) functions that we present 
in Appendix [5] for completeness. Each of these solu- 
tions is composed of the sum of a homogeneous and an 
inhomogeneous solution. The former introduces integra- 
tion constants chosen by requiring (i) that the metric be 
asymptotically flat at spatial infinity, e.g. Hoeo — > as 
r — > 00, and (ii) that the mass and spin angular momen- 
tum associated with the new solution is given by M and 
Ma, as measured by an observer at spatial infinity. 

The line element to 0(a' 2 x' 2 ) is then ds 2 = ds 2 K + 
S(ds 2 ) cs , where 

S(ds 2 ) cs = 2g%dtd4> + 9? t s dt 2 + g^dr 2 + gTede 2 + g^dcf 

(40) 
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with 



ef,' 
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1792 
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4474 M 2 2060 M 3 1500 M 4 
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4221 r 2 469 r 3 469 r 4 

M 2 2220 M 3 1512 M 4 
194— 



2140 M 5 9256 M 6 
+ 



201 r 5 



201 r« 



5376 M 7 
67 r 7 



7 



+ 0(a> 2 X > 4 ) 
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(3 cos 2 6> - 1) 



(41) 



(42) 



(3 cos 2 6 - 1) 



(43) 
(44) 



and g££ given in Eq. (|23|) . Wc have checked explic- 
itly that the solution above satisfies the field equations 
[Eq. (|3"0"|) ] to C(a' 2 x' 2 ) with symbolic manipulation soft- 
ware. 

The choice of homogeneous integration constants de- 
pend on how we choose to define the mass M and the 
reduced spin angular momentum a. The most natural 
choice is to define these quantities as measured by an ob- 
server at infinity, which then leads to the metric displayed 
above. With these definitions, the angular velocity and 
area of the event horizon arc modified to 



9tt 



709 , 
1 - 3^rC , (45) 
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A H , K 1 



915 
28672 



Cx' 



(46) 



where r H is the location of the horizon, which will be dis- 
cussed in the next section (Eq. (|55|l ). and where O h .k = 



a /( r H.K + a 2 ) and A H , K = 16ttM 2 (1 - X 2 /4) + 0( X 4 ), 
the horizon radius for the Kerr metric: r HK = 



One can physically interpret this result 



with 

M + 7M 2 

by thinking of the BH metric found here as representing 
a BH surrounded by a scalar field "cloud" with a dipolc 
density structure. As such, the modified BH horizon area 
and its angular velocity will be modified from what one 
would expect for a Kerr BH in vacuum GR due to the 
presence of the scalar field. 

One could of course insist on defining the mass and 
reduced angular momentum such that the horizon struc- 
ture of the modified BH remains identical to that of the 



Kerr metric. This would require a renormalization of M 
and a via 



- J^CX 2 ) , 



a = a I 1 



57344 
709 



7168 



C , 



(47) 
(48) 



which then leads to = £7 h .k and A H = A H K where a 
quantity with a tilde means that M and a appearing in 
that quantity are replaced by M and a. One can think 
of M and a as a "bare" BH mass and reduced angular 
momentum. The asymptotically defined M and a are 
larger than the bare ones due to the cloud of scalar field 
outside the BH. 

Another quasi-local measure of mass M Komar and re- 
duced spin angular momentum o Komar can be defined by 
the Komar integrals [44|. For the modified BH metric in 
Eq. (gDD one finds 



M Kon 



1 



i ' n[ a rp]^/geeg<t><t>) \ r =r H d9 



M 1- 



1727 
14336 



Cx' 



(49) 



i r w 

— I ((/} a ^n[ a r 0]y /geeg^) \ r =r H d0 



a 1 



29 
128* 



(50) 



where t a d a = —d/dt and 4> a d a = d/dcj> arc timclikc 
and spacelike Killing vectors, while n a and r a are the 
unit covariant vectors normal to the t = const, and r = 
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const, hypersurfaces, respectively. Notice that the above 
quantities are calculated at r = r H , which are different 
from the so-called Komar mass and angular momentum 
which is defined over a 2-sphere at spatial infinity [H| . As 
before, the above quasilocal quantities are smaller than 
those defined at spatial infinity due to the presence of the 
scalar field cloud close to the modified BH. One could rc- 
cxprcss the metric in terms of (M Komar , a Komar ) or (A/, a), 
but we choose not to do so as these quantities are not the 
masses and angular momenta that an observer at infinity 
would measure. 

One might be worried that the solution presented in 
Eq. (|40j) diverges at the unperturbed Schwarzschild hori- 
zon. This is simply a spurious divergence that arises due 
to the slow-rotation expansion, even for the Kerr metric. 
For example, if one takes the (r, r) component of the Kerr 
metric in BL coordinates and expands it in x *C 1, one 
finds terms to leading order that scale as (r — 2M) _1 , 
which diverges at the Schwarzschild horizon. The unex- 
panded Kerr metric, however, can have a horizon located 
inside 2M for prograde spins. One can then be faced with 
the unpleasant situation of the slow-rotation expansion 
of g rr diverging outside the true event horizon due to the 
slow-rotation approximation. 

For practical reasons, one might wish to eliminate this 
feature through resummation. By the latter, we mean a 
modification of certain terms in the metric that naively 
diverge at the Schwarzschild or unperturbed Kerr hori- 
zon, such that 

1. when the resummed metric is expanded in x ^ 1; 
it becomes identical to the old metric to a given 
order in X , i.e. E[g~] = E[g^]. 

2. all components of the resummed metric <?^ c „ sum re- 
main finite everywhere outside the dynamical CS 
modified horizon. 

where g^ c „ sum is the resummed metric, g^ v is the metric 
of Eq. (|40|) and the E[-] operator stands for expansion in 
X< 1. 

In principle, there is an infinite number of ways in 
which one can resum the metric. One way is to re- 
place A — > A' in g rriK [i.e. in the (r, r) component of 
the Kerr metric in Eq. (TTgj) ] and f(r) —> f'(r) in S(ds 2 s ) 
of Eq. (|4T))) , where we have defined 



A 



/» = 1 



915 
" 14336 



CMV: 



(51) 

(52) 



To retain the asymptotic behavior in the x <SC 1 limit, one 
then needs to add the following counterterm in the (r, r) 
component of the metric (induced by the A' modification 
to the Kerr metric): 



2,,2 



915 C,M X 
14336 /'(r) 2 



With these changes, the resummed metric is 

5™ = <W [A'] + [/'] + 5g rr S;s: , (54) 
will satisfy the conditions enumerated above. 

IV. PROPERTIES OF THE SOLUTION 

In this section, we discuss various properties of the 
new solution. We begin by finding the corrected loca- 
tions of the horizon and ergosphere. We continue with 
a calculation of the corrected quadrupole moment of the 
spacetime. We then determine the Pctrov type of the 
new solution. 

A. Singularity, Horizon, and Ergosphere 

The spacetime in Eq. (|40[) contains a true singularity 
at r = 0. This can be verified by computing the Kretch- 
mann invariant -R Aa /p fT -ff' il/pcr , which diverges at r = 0. 
Indeed, this quantity is identical to that found in to 
0(a' 2 x')- We do not present the 0(a' 2 x' 2 ) term here, as 
this cannot cure the r = divergence. 

The location of the event horizon can be found by solv- 
ing the equation gttg<p<t> — g\<k — for r [3]. We find 



915 
28672 



(M X 2 + 0(a u x' A ) 



(55) 



The horizon radius decreases relative to the Kerr horizon 
radius, but of course, the shift of the horizon location 
depends on how one renormalizes the mass and spin and 
also on the choice of radial coordinate. 

The location of the ergosphere can be found by solving 
the equation g tt = for r. We find 



915 

r-cr g o,K - 

+0(a'Y 3 ) 



CM X A 1 + 



2836 
"915" 



■ sin 



(56) 



with the ergosphere of the Kerr solution given by 
r C rgo,K = M + VM 2 - a 2 cos 2 6 fn r H>K [l + (x/4) sin 2 9]. 



B. Lorentz Signature 

We show the Lorcntzian signature of the metric is pre- 
served outside the horizon, provided the coupling con- 
stant is small. Otherwise, our perturbative construction 
of solution will not be justified well. By denoting the 
determinant of the metric component as g and the one 
for the Kerr as g K = —r 2 sin 2 9(r 2 + a 2 cos 2 9) + 0(x' 3 ), 



(53) g/g K is given by 
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g 4211 M 3 2 / 4151/ 14008 M 2 3200 M 3 14780 M 4 8108 M 5 16128 M 6 ^ 

^~ ~ + 6272 r3" X ^ + 201 V + 4221 ~ + 201 T 3 1407 ~ ~ 201 7^ 67" 
201 M 3 , / 1420 M 19888 Af 2 6350 M 3 40100 A/ 4 8524 M 5 16128 M 6 \ 

- — C^Y 2 1 + — 1 7- + — H — r- + — = : tt- 



C r 3 X '( 1 + 603 r + 4221 r 2 + 201 r 3 + 1407 r 4 + 201 r 5 67 r 6 J + ^("'V 3 ) ■ ( 57 ) 



Notice that g/g K does not diverge at r = 2M, just like 
the resummed metric. This is because the determinant of 
the metric is given by g = g rr gee{gttg<p<p - 3^), and while 
9rr oc f(r)~ 2 , the quantity (g tt g^ - 5 t 2 ) oc A 2 - f(r) 2 
and thus g is finite at r = 2M. Since the correction terms 
fall off rapidly as r — > oo, it is important to look at the 
signature of g/g K at the horizon, which is given by 



g 74849 , , / 27901 2 \ , n , ,, ... 

— = 1 Cx 1 + cos 2 9 + 0(a' 2 Y ) ■ 

g K 401408 ^ A V 74849 J y * ' 

(58) 

The correction terms above are negative for any 9. The 
magnitude of the correction to g/g K becomes the largest 
at the poles and at the equatorial plane, respectively. 
One can see that within the small-coupling and slow- 
rotation regime, the signature flip does not take place. 



C. Closed Timelike Curves 

The new BH solution contains no closed timelikc curves 
(CTCs) outside the horizon. If they existed, these curves 
could be found by solving for the region where g^ > 0. 
The explicit forms of g^ at the horizon is given by 



g<t><t> 



4 sin 2 6M 2 



1 



:X cos 



12283 
100352 



Cx 2 



54483 
24566 



cos 2 9 



0{a' 2 x' A ) 



(59) 



where we note that g^ vanishes at the poles. The cor- 
rection terms are positive in the polar region and nega- 
tive in the equatorial region. Eq. (|5"!)|) clearly shows that 
small perturbation due to CS coupling does not change 
the causal structure of spacctimc. 



D. Multipolar Structure 

Since M and a are asymptotic quantities, the first non- 
vanishing correction to the spacetime's multipolar struc- 
ture on g^v appears in the mass quadrupole moment. 
Following Thorne [45[ , we can read off the multipole mo- 
ments by transforming the metric from BL-type coor- 
dinates to so-called asymptotically Cartesian and mass 
centered (ACMC) coordinates (the coordinates where 
the multipole moments are defined in a spacetime region 



asymptotically far from the source). In order to deter- 
mine the quadrupole moment, we need to transform to 
ACMC coordinates such that g u and g^ at 0(r~ 2 ) do not 
contain any angle dependence. In these coordinates, the 
metric component gu for a stationary and axisymmetric 
spacetime can be expressed as 

gtt = + — + [720^20 + {I = P 0lc)]+0 

r 2 r 6 

' (60) 

Here, Y20 is the (£,m) = (2,0) spherical harmonic, and 
I20 corresponds to (m = 0) quadrupole moment. 

Let us first extract the quadrupole moment of a Kerr 
BH. By choosing the flat- spacetime normalized basis 

e = d t , e r = d r , e e =r~ 1 d e , = (r sinfl)" 1 ^ , 

(61) 

associated with BL coordinates, the Kerr metric can be 
re-expressed as 

2M 2Ma 2 cos 2 9 m ( 1 \ 
9u,k = -1 + s + O U , (62) 



9t<p,K 



9rr,K 



gee,K 



g<t><t>,K 



2Masm9 2 Ala 3 sin 9 cos 2 1 

o "T" A 



2M AM 2 - a 2 sin 2 ' 



r3 _ oj,f„2/o _ „„2 , 



8M 3 - 2 A/a 2 (2 - cos 2 9) ( 1 
^ + 



a 2 2Ma 2 sin 2 i 



O Mr 



(63) 

(64) 

(65) 
(66) 



In order to read off the quadrupole moment, we need to 
transform to appropriate ACMC coordinates {if , r', 9', </>') 
such that —a 2 sin 9/r 2 in g rr ^ K and 3ss,k disappears. 
This can be realized by the transformation 



t = t' . 



2r' 



2? 



.12 



(67) 
(68) 

(69) 
(70) 



In this ACMC coordinates, g' tt K becomes 

2M 3A/a 2 cos 2 9' 
r 



' 1 , 

^,k = -1 + -T 33 +°bs)- ( 71 ) 
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Therefore, by comparing this with Eq. 
off the quadrupole moment as (45j 



60l) , one can read 



Here, the quadratic and cubic Weyl quantities / and J 
are defined by [48j 



'20.K 



7T .-, 

—Ma 



(72) 



Let us now follow the above procedure to determine 
the quadrupole moment of the new BH solution. Since 
the correction in the metric is already to G(a' 2 x' 2 ), it 
is not affected by the above coordinate transformation. 
The quadrupole moment in the new solution can then be 
read off as 



ho = I: 



20. K 



1 - 



201 
1792' 



(73) 



Notice that this correction vanishes to linear order in the 
spin; the linear-order in spin terms correct the multipo- 
lar structure of the spacetime at much higher multipolc 
order 0. 

Geroch and Hansen (4(| |4?J proposed a slightly dif- 
ferent definition of multipole moments, which for exam- 
ple leads to a /20,k that differs from Eq. ([72"]) by a fac- 
tor 8yV/15. However, this difference is just a matter 
of convention. One should realize, of course, that the 
quadrupole moment itself is not a directly observable 
quantity. Modifications to the BH multipolar structure, 
however, do imprint on the motion of massive and mass- 
less bodies. Corrections to the gravitational radiation 
induced by this modified motion is indeed observable. 



1 — 2 "^7*^ 



(77) 



J 



with 



3* 2 - 4*1*3 + *4*0 , 
-*2 + 2*i* 3 *2 + *0*4*2 - *4*i - *0*3 



1 / i 



(78) 



(79) 



If Eq. (|76j) is not satisfied, the spacetime is of Petrov type 
I. The Kerr BH and the slowly-rotating BH in dynamical 
CS gravity to 0(a' 2 x') is known to be of Petrov type D, 
which means that Eq. (f75"j) has double degenerate roots. 
In type D spacetimes, not only Eq. (|T6[) holds but there 
are additional conditions that need to be satisfied: 



K = 0, N — 9L = . 



where K, L and N are defined as 
3*4*3*2 



K 
L 

N 



*1*4 
* 2 *4 

*ir 

*4*0 



2*|, 



*3> 



(80) 



(81) 
(82) 



3L 2 



4*4*1*3 



6*4*2*3-3*3. (83) 



E. Petrov Type 

A generic spacetime can be classified into Petrov types 
by counting the number of distinct principal null direc- 
tions (PNDs) W of the Weyl tensor C^ P<J [HEl, where 
fc M satisfies 



k k p k^ T C l j^ l jp[ (T k- x .] — . 



(74) 



This is equivalent to finding the number of distinct PNDs 
i M that make one of the Weyl scalars *o = 0, which 
reduces to counting the number of distinct roots of the 
following equation for b 481 : 



*o + 46*i 



6fe 2 * 2 



46 3 *3 



6 4 * 4 



0. 



(75) 



Here, *o,...*4 are five complex Weyl scalars in an arbi- 
trary tetrad with the restriction *4 7^ 0. 

If Eq (|75|) contains at least one degenerate root, the 
spacetime is said to be algebraically special and the fol- 
lowing relation holds: 



I = 27 J 2 



(76) 



Equation (|76[) determines whether a spacetime is al- 
gebraically special, but when the spacetime is an ap- 
proximate solution, to what order in perturbation theory 
should this equation be calculated? Let us first concen- 
trate on the BH metric in dynamical CS gravity with only 
the odd-parity terms of 0(a' 2 x')- One can construct a 
null tetrad that is a deformation away from the Kerr prin- 
cipal null tetrad, such that * 2 = 0(1), while *i and * 3 
are of 0(a' 2 x') ■ The remaining Newman- Penrose scalars, 
*o and *4, would vanish to this order. One then sees 
that the first term in Eqs. §ff$ and (JTHJ) is of 0(1), the 
second is of C(a' 4 x' 2 ) and the others vanish to this or- 
der. Therefore, Eq. ([75]) is trivially satisfied to C(a' 2 x'), 
while the first non-trivial dynamical CS corrections en- 
ters at 0{a' 4 x' 2 )- Similarly, if one were studying the 
new BH metric found in this paper, i.e. including terms 
of 0(a' 2 x' 2 ), then one would have to consider Eq. ([75)) 
to 0(a' 4 x' 4 )- Obviously, if Eq. fl7j>) is not satisfied at 
0(a' 4 x' 2 ), then one does not need to consider the higher 
order terms. Notice also that the I and J quantities are 
invariant, and one could have chosen another tetrad, but 
the arguments presented above would still hold. 

Now that the order to which terms must be expanded 
is clear, let us focus again on the BH metric in dynamical 
CS gravity with only the odd-parity terms of 0(a' 2 x')- 
Sopuerta and Yunes claimed that this metric is of 
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Petrov type D. We have verified this claim as follows. 
First, we showed that Eq. ([75)1 is satisfied to 0{a' A x' 2 )- 
Then, we showed that the relations in Eq. (|80|) are also 
satisfied to 0{a' A x' 2 )- Although this is sufficient to claim 
that the metric to this order is of Petrov type D, we also 
verified explicitly that Eq. (|75]l has double degenerate 
roots. This then implies that one can rotate the null 
tetrad to a principal one (of the dynamical CS metric), 
where ^ is non-vanishing and contains Kerr terms [of 
0(a' ) but with spin corrections], as well as terms of 
0(a' 2 x')- All other Newman-Penrose scalars vanish at 
this order, i.e. they are at least of 0(a' 2 x' 2 )- 

Let us now focus on the new BH metric in dynamical 
CS gravity, which includes terms of 0(a' 2 x' 2 )- Picking 
the principal null tetrad of the BH metric in dynamical 
CS gravity with only the odd-parity terms included, one 
can show that Eq. ([75)1 is not satisfied to (D(a' A x' A )- One 
might worry that to make this statement precise, one 
would have to account for terms of 0(a' i x') m the grav- 
itomagnctic sector of the dynamical CS metric. These 
terms, however, would modify Eq. ([75]) at 0(a' 6 x' 2 ), and 
thus they can be neglected. Since Eq. ((75)) is not satis- 
fied to 0(a' 4 x' 4 ), the new metric presented in this paper 
breaks symmetries that the odd-parity BH metric used 
to have. This suggests the exact BH solution should be 
of Petrov type I. 



V. GEODESIC MOTION AND SEPARABILITY, 
BINDING ENERGY AND KEPLER'S LAW 

In this section, we discuss the separability of the 
geodesic equations in the modified metric and find the 
binding energy, Kepler's Law and the innermost stable 
circular orbit (ISCO). 



A. Geodesic Motion of Test-Particles and 
Separability 

Consider the motion of non-spinning test particles in 
the new BH solution. We concentrate here on non- 
spinning objects, as otherwise we would have to intro- 
duce an additional scalar dipole-dipole interaction [32| . 
which will be investigated elsewhere [50| . 

One of the most interesting properties of the Kerr 
metric is that the geodesic equations are Liouville inte- 
grate [51(. This leads to the existence of four constants 
of motion or invariants (quantities that Poisson commute 
with the Hamiltonian) : the mass, energy, angular mo- 
mentum and the Carter constant. The existence of this 
last quantity, found by Carter [52 |. as well as the use of 
the proper coordinate system, is crucial in showing that 
the Hamilton-Jacobi equations are separable. In GR, this 
is related to the Kerr solution being of Petrov type D [53| , 
i.e. its associated Weyl tensor possesses double degener- 
ate principal null directions. As shown in Sec. IIVE1 the 
new solution derived in this paper is of Petrov type I 



and there is no guarantee that it possesses a Carter-like 
constant. 

The Carter constant is associated with the existence 
of a second-rank Killing tensor £ Q( 3 , which in GR and in 
BL coordinates is defined by 



tap = Ak( a lp) + r 2 g aj3 



(84) 



with two null vectors k a and l a . The odd-parity BH 
solution in dynamical CS gravity of flol ] does possess a 
Killing tensor at 0(a' 2 x'), where the null vectors k a 
and l a are given by 



k a d a 
l a d a 



A 



-d t + d r +(^-6 g ™)dt, (85) 



(86) 



and 



f cs 5 v M 6 / 12 M 27M 2 \ 



Moreover, one can show that the two null vectors k a and 
l a are also principal null directions of the spacetime to 
0{a' 2 X ') m 

Let us now study whether a non-trivial second-order 
Killing tensor continues to exist to 0(a' 2 x' 2 ), i-e. we 
look for a correction <5£ Q( 3 to the Killing tensor 



top = U + a'V^a/J + 0(a' 2 X ' 6 ) , 



that satisfies the Killing equation V( 7 £ a/ 3) =0. If a con- 
served quantity contains both even and odd parts un- 
der the simultaneous reflection t — > — t and <fi — > —<f>, 
they should be separately conserved. The new BH met- 
ric is symmetric under this simultaneous reflection, and 
hence any geodesic remains geodesic under this trans- 
formation. This means that if we consider a quantity 
£, a @u a uP (here, u a is a four velocity vector), the only 
non-vanishing components allowed are those even in re- 
flection, i.e. (t, t), (t,4>), {(f>, (j>), (r,r), (r,6) and (0,9). 
Without loss of generality, these six components can be 
parametrized through six free functions, A(r, 0), B(r, 9), 
C(r,8), D(r,9), E(r,8) and 5£ r g(r,9), through the fol- 
lowing ansatz: 



5f afl = A{r,e)t a t p +B{r,6)t a <i)p 
+D(r,6)g a p + E(r,6)£ a p 



■C(r,0)<f> a <l>fi 



(89) 



where 8 a d a = d/d6. 

From the symmetry arguments described above, the 
Killing equations contain only 10 independent compo- 
nents. The five functions T = (A, B, C, D, E) appear 
only in the form d r IF or dgJ 7 . Thus, we can solve the 
10 Killing equations for the 10 functions d r T and 8qF 
in terms of S^ r g(r,9). We have found, however, that the 
consistency relation dg(d r B) — d r (dgB) does not hold 
for any <5£ Q/ 3 to 0(a' 2 x' 2 )- Since our ansatz is sufficiently 
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generic, this strongly indicates that there does not exist 
any non-trivial second-rank Killing tensor in the new BH 
solution. 

By ap plying Theorem 1 in Bcncnti and Francav- 
iglia [41| . the non-existence of a non-trivial second-rank 
Killing tensor is enough to claim that the 4 dimensional 
manifold does not admit a separability structure. This 
means that there does not exist any spacetime coordinate 
transformation that leads to the Hamilton- Jacobi equa- 
tion being separable. Since this is an important point, we 
have verified it in two additional ways: (i) by performing 
a Levi-Civita test [42|, and (ii) by trying to map the 
new solution to the most generic spacetime that admits 
a separability structure [41j [see Appendix [B] for more 
details]. In all cases, it is clear that the new metric does 
not admit such a structure. 

Up until now, we showed in various ways that a non- 
trivial second-rank Killing tensor that is a perturbation 
of the Killing tensor found in |4(j does not exist. One 
might wonder whether there is a completely new second- 
rank Killing tensor that is not a perturbation of that 
in [Z(|. If this exists, one must be able to find it by set- 
ting £ a p — in Eq. (|88[) . By imposing the same ansatz 
of Eq. (|89|) . we can again solve for d r JF and dgJ- ' . This 
time, dg(d r B) = d r {dgB) is trivially satisfied. However, 
one finds that the only solution of 8^ r g that satisfies 
dg{d r F) = drldgJ 7 ) is 6£ r g = 0, leading to <5£ Q( 3 = 0. 
This proves that a completely new non-trivial 2nd-rank 
Killing tensor at 0(a' 2 x' 2 ) does not exist. The non- 
existence of this tensor can be proved in a different man- 
ner. If it exists, the leading contribution should start at 
0(a' 2 x' 2 )- However, the only possible form of this lead- 
ing term would be ct' 2 x' 2 ^afj wnere £u/3 * s * ne 2nd-rank 
Killing tensor of the Kerr spacetime. This is because the 
completely new Killing tensor divided by a' 2 x' 2 should 
also be a Killing tensor whose leading term should satisfy 
the GR Killing equations. Since ct' 2 x' 2 (,ai3 does not sat- 
isfy the Killing equations at 0(a' 2 x' 2 ), we conclude that 
the completely new 2nd-rank Killing tensor cannot exist. 
Of course, we cannot rule out the possibility of the exis- 
tence of the completely new Killing tensor once the small 
coupling or slowly rotating approximation is violated. In 
this sense, we have only shown the non-existence of the 
perturbative non-trivial, 2nd-rank Killing tensor in the 
new BH solution. However, we emphasize that, in this 
paper, we only focus on the situation where both of the 
approximations hold. 

Although the geodesic equations are not exactly inte- 
grable, the new solution is sufficiently close to the one 
found by Yunes and Pretorius Q that, except for the 
resonant orbits, the geodesic equations are still approx- 
imately integrable. By the latter, we mean that when 
one orbit-averages, there still exists a Carter-like con- 
stant, i.e. the 0(a /2 x' 2 ) terms that spoil the existence 
of a Killing tensor are odd in ut, where u is any of the 
fundamental frequencies of the motion, and thus, van- 
ish upon orbit-averaging. This can be shown explicitly 



by applying canonical perturbation theory [55[ following 
e.g. Glampedakis and Babak [5(|, as we discuss in Ap- 
pendix [Cj 

The new metric found here cannot be mapped to the 
new bumpy metrics proposed in 57|. This is because 
the latter assumed the existence of a non-trivial second- 
rank Killing tensor, while the solution found here does 
not possess it. We have tried to map the new solution to 
a generic deformed Lewis-Papapetrou spacetime [58| in 
one of the BL-typc coordinates. We found that a naive 
mapping does not seem to work, which implies that a 
further coordinate transformation is probably needed. 



B. Binding Energy, Kepler's Third Law, the 
Location of the ISCO and Curves of Zero Velocity 

From the definitions of the energy E and the (z- 
component of) orbital angular momentum L z , we have 



Eg^ + L z g t 4, 

9t4> 2 - gttg<t><t> 
Eg t(j> + L z g tt 

g 2 tt p - gugu 



(90) 
(91) 



where the overhead dot stands for a derivative with re- 
spect to the affine parameter. By substituting the above 
equations in u a u a = —1, with u a the particle's four- 
velocity, we obtain 



g rr r 2 + g ge 6 2 = V Btt (r, 6; E, L z ) . 



where the effective potential is given by 

?2„,, i opt n . . _i_ r2. 



V, 



_ E g<p<p + 2EL z g t(j> + L 2 z g tt 

it — 9 ^ 



% - gug<p<p 



(92) 



(93) 



For simplicity, we restrict attention to equatorial, cir- 
cular orbits. Then, E and L z can be obtained from 
V c!! = and dV c!! / dr = as 



E = E K + SE , 

L z = L z ■ K + 5L Z 



(94) 
(95) 



Here, E K and L ZjK are the energy and the (z-componcnt 
of) orbital angular momentum for the Kerr back- 
ground [5l| : 

r 3/2 _ 2M r i/2 + aM 1 ' 2 
K = r 3 / 4 (r 3 / 2 - 3MrV2 + 2aM 1 /2)i/2 ' ( > 

M l/2( r 2 _ 2aM l/2 r l/2 +Q2) 
Z ' K = r 3/4( r 3/2 _ 3Mr l/2 + 2aM 1 /2)l/2 ' ( ) 

where we have defined <f> to be positive in the direction of 
prograde orbits. This implies that negative a corresponds 
to retrograde orbits. The CS corrections are 
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5 M 11 / 2 f 33 M 183 M 2 603 A/ 3 81 A/ 4 \ 

= 4^ x r 3( r „ 3M )5/2 ^ _ TIT" ~ lZoT 2 " ~ To"" + T~) 

201 2 M 3 / M 59315 M 2 38954 Af 3 289564 A/ 4 188420 M 5 

~7168^ X r 1 /2( r _3Af)5/2 I + ~ ~ 42 2l r 2 + 4221 ~ + 4221 T 4 " + 1407 r 5 

566500 A'/ 6 27360 Af 7 61584 A/ 8 96768 M 9 



1407 r 6 67 r 7 67 r 8 67 r 9 



+ 0(a' 2 X ' 3 ), (98) 



15 A/ 5 / M 2M 2 A/ 3 108 M 4 

z = "8" Cx r 3/2( r _ 3M )5/2 ^ _ —-572"- 6 7T + — 7F 

603 A 2 rM 5 / 2 AM 54833 M 2 , 110798 A/ 3 . 15100 M 4 . 369428 M £ 



;CX 2 7 ^tt^ 1-77 — 



7168^ A (r - 3Af) 5 / 2 ^ 3 r 12663 r 2 12663 r 3 12663 r 4 4221 

74092 A/ 6 32768 M 7 40688 A/ 8 32256 A/ 9 ^ { ^,2^,3 
603 r e + 469 T 7 67 r 8 " + ~~67 r 



+ 0(a'Y 3 )- (99) 



When we expand E and L z in powers of M/r, the 
leading-order correction to the binding energy E b = E— 1 
and L z are 

/ 201 , , A/ 2 \ / /2 Af 4 

(100) 

Relative to the leading-order Kerr (or Kepler) terms, the 
corrections are proportional to (M/r) 2 which are of 2PN 
orders. As before, the corrections in E and L z would 
change if one used a different renormalization of the mass 
and spin, such as M and a, but these quantities are not 
observable at spatial infinity. 

We can also derive the correction to Kepler's Third 
Law by calculating the orbital angular frequency of a 
test particle w = L z /r 2 to find 



u 2 = u>l I 1 



where w 2 = A/(r 3 / 2 + aA/ 1 / 2 )" 2 [Hj]. 

However, the expressions for E, L z and uj above are not 
gauge invariant. The gauge invariant relation between E 
and u> can be obtained by expanding Eqs. (jlOOp and (|102j) 
to 2PN order and eliminating M/r. The final result is 



oj(E) 



2V2 

~aT 

891 

+ *r 

+0 [\E b 



E h \ 3 ' 2 



64 



\E b \-AV2\E b 



67 



8/2 



297^ 2772 
±0(a' 2 X ' 3 ) 



C% 2 )\E b \< 



(103) 



and its inverse is 



E(lu) = 1--(Mlj) 2 / 3 + -(Mu J ) 4 / 3 -^(Mlo) 5 / 3 
2 8 3 



27 
~16 

-O 



1+ 27 

(Afw) 7 / 3 



V 



1 67 
2016 

+ CVY 3 ), 



CX 2 (Ma;) 2 



(104) 



To O(a'°x /0 ), this agrees with the standard PN re- 
lation shown in j59| . 

Let us now derive the correction to the location of 
ISCO. Substituting Eqs. (JMJ and (JS5J) in Eq. ^93]) , and 
then solving the equation d 2 V<, s /dr 2 = for r, we obtain 



77^6 



5184 



9497219 



219469824 



(M X 2 



(105) 



^ISCO — 7~ISCO,K 

+0(«' 2 X' 3 ): 

where the Kerr ISCO radius is given by [5"i| 

r ISCO , K = A/ {3 + Z 2 - [(3 - Zi)(3 + Zi + 2Z 2 )] 1/2 } 

(106) 

with 

^1 = i + (i-x 2 ) 1/3 [(i + x) 1/3 + (i-x) 1/3 ], 

(107) 

Z 2 = (3 X 2 + Zj) 1 / 2 . (108) 

The CS correction at linear order in \ agrees with that 
found in while the 0(a' 2 x' 2 ) term is new. The radial 
location of the ISCO, however, is not gauge invariant. A 
gauge invariant quantity can be obtained by calculating 
the angular orbital frequency uj 1sco at ISCO, which is 

77 



w 1: 



CJr, 



X 



124416" M 

,2 



(109) 



31603654656" Af 



14 



where w ISCO 



ncic u^ ISCO — ™ ' V' isccxk t x lyl I 

For completeness, let us also compute the correction 
to the radiative efficiency rj, which is defined by 



T)=l — E(r ISCO ) 



(111) 



This quantity corresponds to the maximum fraction of 
energy being radiated when a test particle accretes into 
a central BH. For Schwarzschild and extremal Kerr BHs, 
rj ~ 0.06 and r\ ~ 0.42, respectively. The radiative effi- 
ciency for the new solution is then 



31352832 



48771072 



Notice that there are both linear in x an d quadratic in 
X corrections. 

Finally, let us consider curves of zero velocity 
(CZV) [60, [6l| in the r — 9 plane, i.e. curves with V ef! = 0. 
Since the left-hand side of Eq. is always positive, 
bound orbits are allowed if and only if V etf > 0. Fig- 
ure [T] shows CZVs for the Kerr solution while Fig. [5] 
shows the one for the new solution. The enclosed regions 
represent the regions where V eH > and the solid lines 
at r/M = 1.95 correspond to the location of the event 
horizon for the particular case considered in the figures, 
i.e. E = 0.95, L z = 3M, X = 0.3 and in the CS case, 
C = 0.1. When we draw these figures, we first expand 
the metric g M „ in a and calculate V eff . We do not further 
expand this V att in a, since if we do this, V ctt is propor- 
tional to a negative power of a, which would render the 
expansion invalid near the horizon. 

Two allowed regions are clearly visible in these figures: 
an outer region and an inner region, for both the GR and 
CS cases. The outer regions are similar in GR and CS 
theory, although due to the scale of the figure the differ- 
ences look small. As shown in [4(| and recently in plj . 
orbits in the outer region are still distinguishable with 
gravitational wave observations. On the other hand, the 
structure of the inner regions change drastically, as ex- 
pected since the CS correction modifies the strong field 
regime. These inner enclosed regions arc inside the hori- 
zon, however, and thus they cannot be probed with grav- 
itational waves, at least for slowly rotating BHs. Also, 
we cannot trust the perturbative solution there. 



VI. CONCLUSIONS AND DISCUSSIONS 

We have found a stationary, axisymmetric BH solution 
in the small-coupling and slow-rotation approximations 
at linear order in the coupling constant but at next-to- 
leading order in the spin. This solution does not satisfy 
the vacuum Einstein equations but the modified field 
equations. We used a novel technique to find this so- 
lution, based on Schwarzschild BH perturbation theory. 
That is, we decomposed the metric perturbation and the 
source terms (that come from modifications to GR) in 
tensor spherical harmonics, reducing the field equations 



to a set of coupled, ordinary differential equations that 
are much simpler to solve. We found that corrections 
at quadratic order in the spin appear in the even parity 
sector of the metric. 

The method presented here could be used to find so- 
lutions both to higher order in x (or a/M) and to higher- 
order in £ defined in Eq. ([3]). The dynamical CS action, 
however, is a linear-order-in-£ truncation of a more fun- 
damental theory, and thus, it is valid only to linear order 
in the coupling constant. If one were to carry out this cal- 
culation to O(o. x')> one would find a modification only 
in the gravitomagnetic sector, which is easy to compute. 

A nontrivial property of the new solution is that, al- 
though it is of Pctrov type D to 0(a' 2 x'), it is of Petrov 
type I to 0(a /2 x' 2 )- This is different from the Kerr met- 
ric, which is of Petrov type D to all orders in x' ■ The 
new metric does not possess a second-order Killing ten- 
sor or a Carter- like constant. This then implies that 
there does not exist a spacetime coordinate transfor- 
mation that leads to Hamilton- Jacobi equations being 
separable, which also implies that geodesic motion is, 
in all likelihood, chaotic when corrections of 0(a' 2 x' 2 ) 
are included in the metric. However, although there 
is no exact Carter-like constant, we have also showed 
that the geodesic equations are still separable after orbit- 
averaging, except for the resonant orbits, by applying 
canonical perturbation theory (55l. [56l lo2j . In some sense, 
then, it might be possible to recover geodesic regular- 
ity on average, although it is not clear what this means 
precisely. Possible future work could concentrate on 
studying whether geodesies in this background are truly 
chaotic, and if so, whether such chaos manifests itself 
outside the event horizon. Moreover, one could also in- 
vestigate how large the shifts in the orbital frequencies 
of geodesic test-particles are and discuss observational 
prospects of probing such a spacetime (see Vigeland and 
Hughes [62[ for similar work on bumpy spacetimes, as 
well as Appendix [C]). 

Some insight might be gained by comparing geodesic 
orbits in this new background to those in the Manko- 
Novikov (MN) spacetime (63|. Gair et al. [6(| inves- 
tigated geodesic orbits in the MN background through 
CZVs in cylindrical coordinates (t,p,z,<fi). They found 
that when the quadrupole moment deviates from Kerr, 
chaotic islands arise in the inner region of the p — z plane. 
Strictly speaking, the Hamilton- Jacobi equations are not 
separable in the MN metric, but in the outer region, it 
seems that there exists a nearly invariant quantity that 
corresponds to a Carter-like constant. Such a result is 
related to the Kolmogorov, Arnold and Moscr (KAM) 
theorem (5|| which states that when a separable Hamil- 
tonian system is weakly perturbed, the perturbed motion 
within the phase space remains mostly in the neighbor- 
hood of the invariant tori (see related work by Aposto- 
latos et al. [63 - [6r| ). This suggests that except for certain 
resonant orbits, there should exists a fourth constant of 
motion in the phase space where motion is non-chaotic. 
This constant of motion, however, is not related to the 
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£ = £ = 




r/M r/M 



FIG. 1: Curves of zero velocity, V e g = for the Kerr metric. The enclosed regions show the allowed orbit region Vcu > for 
E — 0.95, L, — 3M and \ = 0-3- The left panel corresponds to the inner region while the right panel corresponds to the outer 
region. The thick solid lines at r/M = 1.955 correspond to the location of the event horizon when \ = 0.3. 



£ = 0.1 £ = 0.1 £ = 0.1 




r/M r/M r/M 



FIG. 2: Curves of zero velocity, V„a = for the new BH metric with £ = 0.1. As in Fig. [T] the enclosed regions show the 
allowed orbit region V B tt > for the same parameters. The middle panel corresponds to the inner region while the right panel 
corresponds to the outer region. The left panel zooms into the region around 6 ~ (3/4)7r. Similar structure appears at around 
9 ~ (1/4)tt. The shaded areas are the allowed regions where test-particle orbits exist. The thick solid lines at r/M — 1.955 
correspond to the location of the event horizon when % — 0-3. 



symmetries of the spacetime anymore, and since it does 
not exist for resonant orbits, one cannot expect to find 
a global Killing tensor. In Fig. [TJ we have shown that 
there are no additional CZV islands produced in the new 
spacetime compared to Kerr outside the event horizon. 
Although the structure of the inner region changes drasti- 
cally, this occurs inside the horizon. Whether chaotic or- 
bit exists in the new spacetime outside the horizon needs 
to be verified numerically 

In addition to the global structure of the spacetime, 
we investigated other properties in detail. We first found 
that the horizon radius and the location of the ergosphere 
are modified from the GR expectation, although they are 
not modified up to linear-order in the spin. We also com- 
puted the 2PN conservative corrections to the binding 
energy and Kepler's Third Law, in terms of the asymp- 
totic mass M and spin a. This, in addition to the dis- 
sipativc corrections found in [32j, will allow for the con- 
sistent calculation of the gravitational waveform to 2PN 



order, after the dipole-dipole scalar force is calculated. 
This waveform can then be mapp ed to the parametcrized- 
post-Einstcinian framework [671 ] and future gravitational 
wave constraints on £ can be investigated. 

One might wonder whether constraints on £ can be 
placed with observations of the orbital decay of binary 
pulsars. The gravitational fields outside of a BH, how- 
ever, is very different from that outside a neutron star 
in dynamical CS gravity [H, 113, due to the different 
boundary conditions used when solving the modified field 
equations. Thus, the new BH solution found here cannot 
be used to investigate binary pulsar constraints. Instead, 
one would need to obtain NS solutions to quadratic or- 
der in spin, so that the dipole scalar charge that sources 
2PN scalar radiation and the dipole-dipole force can be 
calculated. This would lead to different 2PN conserva- 
tive corrections to the binding energy and Kepler Third 
Law. 

Two low-mass X-ray binary (LMXB) systems contain- 
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ing BHs have yielded a measurement or constraint on the 
orbital decay [(39|, but the BH spins have not yet been 
determined. Since the dipole scalar charge is sourced by 
BH spins, it is currently impossible to place a constraint 
on dynamical CS gravity with these systems. Moreover, 
binary parameters in binary pulsars and LMXBs have 
been determined by assuming that GR is correct. Since 
the conservative part (that determines the binary param- 
eters) appears at the same order as the dissipative part 
(that causes the orbital decay rate), one would have to 
redo the fits to simultaneously determine binary param- 
eters and constrain f . 

Electromagnetic radiation from accretion disks around 
a central BH, such as Sgr A*, can also be a powerful 
tool to test GR [l2|. One could study how these ob- 
servables, e.g. images of BH shadows [7(1 [Zll, contin- 
uum spectrum f72l |73| , quasi-periodic oscillations [74| , Fe 
emission lines [75| . geodetic precessions and strong lens- 
ing (7(|, are modified if the central super- massive black 
hole (SMBH) is described by the new solution found in 
this paper. We expect that if the central BH is spinning 
moderately fast (though not extremely fast since then 
our slow- rotation assumption breaks down), the CS cor- 
rection at quadratic order in spin will be dominant over 
the linear-order, CS corrections. 

Another possible avenue for future work is to relax the 
slow-rotation approximation. Obtaining an arbitrarily 
fast rotating BH solution analytically seems difficult, but 
a numerical solution might be feasible. The results in this 
paper can then be used as a check of the slow-rotation 
limit of such a numerical solution. Once we obtain the 
correction to the quadrupolc moment without applying 
the slow-rotation approximation, we can map this to a 
correction to the gravitational waveform. 

As a final remark, one could also investigate other 
modified gravity theories with the techniques developed 
in this paper. For example, in Ref. [32J, we studied 
Einstcin-Dilaton-Gauss-Bonnet theory, where a static, 
spherically symmetric BH solution [10[ is known, as well 
as a stationary, axisymmetric BH solution at linear or- 
der in spin [ll[ . We have checked that this solution can 
be mapped to Benenti and Francaviglia metric shown in 
Appendix IB 21 suggesting that there exists a non-trivial 
second-rank Killing tensor and a Carter-like constant, 
and hence it admits a separability structure. One could 
then, for example, check whether the latter is of Petrov 
type D. One could also extend this solution to quadratic 
order in spin to see how the properties of the spacetime 
change compared to the linear order corrections. 
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Appendix A: Tensor Harmonics 

In this paper, we used the following tensor spherical 
harmonics to decompose the metric perturbation and the 
source term [H, [HI : 
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where Yeo are the m = spherical harmonics and W^q 
are given by 
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On the other hand, the coefficients of the source after 
a tensor spherical harmonics decomposition are 
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By substituting the above source terms in Eqs. (|3"5|) - (|3"9")l , we obtain a set of ordinary differential equations for 
(ffooo, #200, K 00 , H 020 , H 2 2o,K 20 ), which we solved to find 
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r 



This solutions can then be used to reconstruct the metric 
perturbation, presented in the main text. 



Appendix B: Alternative Ways to Prove the 
Non-admittance of a Separability Structure 

In Sec. IV Al we have shown that a 2nd-rank Killing 
tensor does not exist for the new BH solution by directly 
attempting to solve the Killing equations with a general 



ansatz for the Killing tensor, restricted by the symmetries 
of the spacetime. This implies that the new BH solution 
does not admit a separability structure. In this appendix, 
we verify this point in two additional ways. 



1. Levi-Civita test 

Whether a Hamiltonian system is separable or not in 
a given coordinate system x a can be determined by the 



18 



so-called Levi-Civita test [42|, |54|. This test states that 
the Hamilton- Jacobi equations are separable if and only 
if the Hamiltonian H = {l/2)g ttv p lJi p v satisfies the 
following equations for each pair of distinct indices a 
and (3: 

dH dH d 2 H dH dH d 2 H 

dx a dx@ dp a pp dp a dpp dx a x@ 
dH dH d 2 H dH dH d 2 H 

= 0. (Bl) 



dp a dx@ dx a pp dx a dpp dp a x fi 

Here, p a is the conjugate momentum of the coordinate 
x a . Note that the repeated indices do not indicate sum- 
mation. For the stationary and axisymmctric spacctimc, 
since the metric only depends on r and 9, the only rel- 
evant combination of (a, /3) in the above equation is 
(a, f3) = (r, 9). We have checked that the slowly rotating 
BH solution in dynamical CS gravity at linear order in 
spin satisfies Eq. (|B1[) . The new solution in BL-like coor- 
dinates does not satisfy the above equations, as we veri- 
fied explicitly with symbolic manipulation software. This 
suggests that there is no conserved quantity correspond- 
ing to the Carter constant that can be constructed from a 
2nd-order Killing tensor, and thus, the Hamilton- Jacobi 
equations are not separable in BL-like coordinates. 

Next, we investigated whether there exist any coordi- 
nates in which the new BH solution satisfies Eq. (|B1I) . 
In particular, wc allowed for diffcomorphisms x^ — ► 



restricted to £^ being of order 0(a' 2 x' 2 ^ 



With this coordinate transformation, g^ iv transforms via 
9jii!/ g'au = 9tiv — 2/i( M;!y ). Wc verified, however, that 
the transformed Hamiltonian does not satisfy Eq. (|B1[) 
with (a,/3) = (r,9). We found that the left- hand- 
side of this equation contains only five terms, which are 
proportional to the combinations {E 2 p r pg), (EL z p r pg), 
(L 2 p r pe), {pipe) and (p r Pg)- If Eq. (jBip is satisfied, the 
coefficient of each term must vanish separately (i.e. there 
are five equations to be satisfied). On the other hand, we 
found that £ p appears in these coefficients only through 
deC ', d r £ e , d r dgt; r and d r de£, 6 . We tried to solve the five 
equations for these four quantities, but could not obtain 
a consistent solution, which then proves that there does 
not exist any coordinate transformation where Eq. (jBip 
is satisfied. 



2. Mapping to the General Metric that admits a 
Separability Structure 

The inverse of the metric components of a spacctimc 
that admits a separability structure can be expressed 
as E3 



where p = acos9. There are four functions of r, Q(r) 
and Cr & ( r )> an d f° ur functions of 9, P(p) and Cp b {p)- The 
Kerr solution can be expressed as 

Q K (r) = A, P K (p) = a 2 sin 2 0, 
1 f{r 2 +a 2 ) 2 a(r 2 + a 2 ) 
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The deviation from the Kerr solution can be 
parametrized as 



2 -(i + *Q), 

(r 2 + p 2 )a 2 sin 2 
1 



r 2 + p 2 



(B4) 
(B5) 



[Q*(% K (l + 6Q + 8(?) 
+PkC"k(1 + W + 0] > (B6) 



9 H = [QkC* k(1 + SQ + 5^) 

7™ p 



where, 



SA = 



A-A K 
A K 



(A 



p j-ab /-a>b\ 



(B7) 



(B8) 



(B9) 



This can be interpreted as the most general bumpy space- 
time that admits a separability structure. 

Let us now try to map the new BH solution to the 
above bumpy metric. First, to 0(a' 2 x'), the only rel- 
evant parameters are SQ, SP, S(* 1 and ^C*^- We find 
that 

SQ = SP = SC = , 

To C(Cx 2 )) we nn d that the new BH solution cannot be 
mapped to the above bumpy spacetime with BL-like co- 
ordinates. We also considered coordinate transforming 
the new BH solution to find a map to the above metric. 
Once again, however, wc found that there does not ex- 
ist any coordinates where the new BH solution can be 
mapped to the above bumpy spacetime. This then im- 
plies that new BH solution does not admit a separability 
structure. 



r = 9 6 



O I o * 



p(p) 



(r 2 + p 2 )a 2 sin 2 9 ' 



ab = Q( r ) r ab ( , , P (P) sab 
I r 2 _|_ p2 V J ~ 7 ,2 j_ „2 



(B2) 



Appendix C: Canonical perturbation theory and 
secular separability of the geodesic equations 

Consider a Hamiltonian system in angle-action vari- 
ables, where we use w a to denote the angle variable and 
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I a = (pt, J r , Js, Jtj>) to denote the conjugate action vari- 
ables, with 



Ji = f dx l pi . 

We can then express the Hamiltonian as 
1 1 



(CI) 



H 



-eh^ Pfl p„ + 0(e 2 



= H {Ip) + eH 1 {w a ,Ip)+O{e 2 ), (C2) 

where g Y p denotes the stationary, axisymmetric CS BH 
metric at 0(a' 2 x') found by Yuncs and Pretorius @ 
and h?" is the difference between the new solution and 
theirs. We have also defined Hq = (\/2)g^pp IJj p ll and 
H\ = {l/2)Wp^p v , where e = O(x') is a bookkeep- 
ing parameter that labels the order of the perturbation. 
Clearly, the first (second) term is the unperturbed (per- 
turbed) Hamiltonian. Note that the unperturbed Hamil- 
tonian depends only on the action variables, while the 
perturbed one depends both on the angle and action vari- 
ables. 

Now, let us seek a canonical transformation from 
(w a ,Ip) to (w a ,I[j) such that the new action Jg is con- 
stant and the new angle variables w a are linear in the 
affine parameter. With this choice, the new Hamiltonian 
would depend only on the new action variables. Such a 
canonical transformation exists if one can find an appro- 
priate generating function, which can be parameterized 
as 15a. l56i 



F(/,y=/I a |*(»Vp) 



(C3) 



where $ is a function of the old angle variables and the 
new action variables, which are assumed to be periodic 
in w a . 

Whether such a canonical transformation exists then 
depends on whether $ exists. This function must satisfy 
the Hamilton- Jacobi equation 



(C4) 



where H is the new Hamiltonian. We expand H as 

H{ip,e)=H {Ip) + eH l {ip)+O{e 2 ), (C5) 
and substitute this into Eq. (|C4|) to yield 

H (I P ) = H (ip), (C6) 
rf-j^+H^w",!?) = H^Ip), (C7) 

where v c q = dHo/dI a represents the unperturbed fre- 
quencies. Next, since $ is assumed to be periodic in w a , 
we Fourier decompose it 



(C8) 



where j is a 3-dimensional vector of the integer in- 
dices. The coefficients Bj can be determined by solv- 
ing Eq. (|C7p. Since the derivative of $ with respect to 
w a does not contain any constant term, the first term 
on the right hand side of Eq. (|C7[) does not depend 
on If; when one orbit-averages. Therefore, upon orbit- 
averaging, Eq. (|C7[) becomes 



(H 1 )(w a ,I )=H 1 (i ), (C9) 
where the orbital averaging is defined as 



(A) = & dw a A , 



(CIO) 



for any quantity A. The oscillatory part of Eq. (|C7I) must 
satisfy 



<9$ 

dw a 



= (Hi) - Hx 



(Cll) 



Similar to Eq. (|C8[) , we decompose the right hand side of 
the above equation as 

{H 1 )-H 1 = Y, Cj(ip)e 2 ^- W . (C12) 



By substituting this equation and Eq. (|C8[) into 
Eq. (jClip . we can determine the coefficients Bj in 



Eq. E8J) as 



(C13) 



Thus, when j ■ i/ ^ 0, the Bj coefficients exist, which 
then implies that $ and the canonical transformation in 
question indeed exist. 

If the generating function exists, the relation between 
the old and the new action variables is 



<9$ 



(C14) 



dw a dw a ' 

Taking the orbital average of this equation, wc obtain 



(C15) 



Thus, there is no secular change in the action variables, 
which means that the conserved quantities (energy, an- 
gular momentum and the Carter-like constant) are not 
modified in a secular sense (56J. Therefore, when we take 
an orbital average, the geodesic equations must remain 
secularly integrable. 

This, however, does not mean that there are no secular 
changes in the angular frequencies u) % . The shifts in these 
frequencies co njug ate to the orbital proper time can be 
calculated via 1621 



Sw* 



m dli 



(C16) 
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where m is the mass of the test particle. This can be 
related to the shifts in the observed angular frequencies 
f2* via 



with 



1 dH 
m dL 



1 dH 



l 5T 



m dL 



V 2 



5T 



1 fl(ffi) 



(C17) 



(C18) 



Resonant orbits must be treated separately. In this 
case, j ■ Vq =0 and the right-hand side of Eq. (|C13|) 
diverges. This, in turn, means a canonical transformation 
does not exist, and thus, a Carter-like constant of motion 
does not exist either at resonance, even when one takes 
the orbit averaging. Since a constant of the motion must 
be a global quantity, this implies that a Carter constant 
cannot globally exist. 
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